Abstract. We report 0.9-1.4 fim spectroscopic observa- 
tions of NGC1068 collected during the commissioning 
phase of the near infrared camera spectrometer (NICS) 
of the Telescopic Nazionale Galileo (TNG). These yielded 
the first extragalactic measurement of [PII] (1.188 fim) 
line emission. In the central 0.75"x2" the [FeII]/[P II] Hne- 
intensity ratio is close to unity, similar to that measured in 
the Orion bar and a factor of ^20 smaller than in super- 
nova remnants. This indicates that most of iron is locked 
into grains and, therefore, argues against shock excitation 
T-H [ being the primary origin of [Fe II] line emission in the cen- 
O tral regions of NGC1068. 

We propose the [Fell] /[PII] ratio as a simple and effec- 
tive tool to study and perhaps resolve the long debated 
I questions related to the origin of [Fe II] line emission and, 
more generally, to constrain the role of shock excitation 
in active galaxies. 
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1. Introduction 

Since the first infrared spectroscopic observations of galax- 
ies and supernova remnants in the 80 's, the emission lines 
of [Fe II] have become a popular and debated issue (Moor- 
wood & Oliva |1988|, Forbes & Ward |1993|, Simps on et 



Mouri et al. |200CD 



al. |1996| , Vei lleux et al. |l997| , Alonso-Herrero et al. |l997 



From the observational point of view, [Fe II] is weak in 
HII regions and planetary nebulae while extremely strong 
in shock-excited filaments of supernova remnants. Since 
relatively bright [Fe II] emission is commonly found in the 
IR spectra of normal and active galaxies, many authors 
have considered the possibility of using this line as shock 
tracer and, even, to count the number of supernova rem- 
nants (e^ Colina 1993 , Vanzi & Rieke 1997 , Engclbracht 
et al. [1998[ ). 

From the theoretical point of view, a low density region 
with normal abundances can become a strong source of 
[Fell] only if the following conditions are satisfied. 

i) Most of iron must be in the gas phase, i.e. dust grains 

must have been destroyed. 
ii) The gas electron temperature must be large enough 
to collisionally excite the upper levels of the lines, in 
practice Z 5000 K. 
in) Most of iron must be singly ionized, i.e. Fe^/Fe must 
be close to unity. 

Given the low ionization potential of Fe"*" and the high ef- 
ficiency of the Fe -I- H° charge-exchange recombination 
reactions, the latter two conditions are equivalent to say- 
ing that bright [Fe II] lines can only be formed in regions 
where hydrogen is partly ionized (e.g. Oliva et al. 1989). 



The most efficient mechanisms for creating extended 
regions of hot, partially ionized gas are shocks and pho- 
toionization by soft X-rays. In both cases the volume emis- 
sion measure of the partially ionized region could easily 
exceed that of the fully ionized gas. The only important 
difference between the two mechanisms is that photoion- 
ization is unable to destroy the toughest iron-based grains 
which are otherwise easily sputtered by shock fronts. 

In practice, therefore, a purely photoionized region 
emitting [Fe II] can be easily recognized by measuring the 
abundance of Fe+ relative to any non-refractory species 
which forms in the same partially ionized region. Finding 
a very low iron abundance would unequivocally imply that 
the gas has not been significantly processed by shocks. 

This could be in principle obtained by comparing the 
infrared [Fell] and optical [01] lines but, in practice, the 
large difference in critical densities, reddening and the 
problems of comparing data taken with different instru- 
ments makes it impossible to obtain a clear-cut conclu- 
sion (see e.g. Mouri et al. |1993l Alonso-Herrero et al. 



1997, Larkin et al. 1998| ). A much more reliable deter- 
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mination of the iron relative abundance should be derived 
from lines close in wavelengths and with similar critical 
densities such as the [P II] and [Fe II] lines which are dis- 
cussed in this Letter. 



2. Observations and results 

The data were collected at the Telescopio Nazionale 
Galileo (TNG) in November 2000 during the commis- 
sioning phase of NICS, the near infrared camera and 
spectrometer expressly designed and built for this tele- 
scope. This instrument is a FOSC-type cryogenic focal 
reducer equipped with two interchangeable cameras feed- 
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Fig. 1. NICS-TNG spectrum of the Seyfert galaxy NGC1068 compared with spectra of template objects from the 
literature (see text Sect. ||). The break at /xm corresponds to the region of bad atmospheric transmission. The 

flux scales and zero levels have been adjusted to facilitate the direct comparison between the spectra. 



ing a Rockwell Hawaii 1024^ array. The camera used for 
the spectroscopic observations has a projected scale of 
0.25'7pixel (Oliva & Gennari [l995| Baffa et al. pOO0| ). 
The spectroscopic modes are achieved by means of a se- 
ries of glass-resin grisms which can be inserted in the 22 
mm collimated beam (Vitali et al. 1997). The spectrum of 
NGC1068 was collected through a sht of 0.75" (=3 pix- 
els) width and using the IJ grism which yields a 0.89- 
1.46 fim spectrum with a dispersion of 5.7 A/pix. The slit 
was oriented N-S (i.e. at PA=:0°) and centered on the 1 
fim continuum peak. The acquisition consisted of a series 
of four 5-minute exposures with the object set at differ- 
ent positions along the slit followed by halogen flats. The 
atmospheric spectral response and the instrumental efS- 
ciency were determined using spectra of the 06.5V star 
HD42088 whose intrinsic spectrum was approximated by 
i^A=l-7-10"^ (A/1. 25)^3.7 gj.g ^^-2 g-i ^i^-i^ 

The spectrum of the central 0.75"x2" region is dis- 
played in Fig. |l] where we also show, for comparison, spec- 
tra of the Orion Bar (Walmsley et al. 2000) and unpub- 
lished spectra of supernova remnants collected in 1992 us- 
ing IRSPEC at the ESO-NTT telescope. The relative hue 
fluxes are summarized in Table |l|. Evident is the differ- 
ence between the very large [Fell] /[PII] ratio measured 
in SNR's and the much smaller values found in the Orion 
Bar and in NGC1068. 



In principle, the emission feature peaking at 1.188 fim 
could be contaminated by [Nill] 1.191 /im. However, this 
line was measured in the Crab nebula at a level of only 



<15% of [Fell] (Rudy et al. |1994| ) and was not detected in 
the supernova remnants listed in Table ^. Moreover, [Ni II] 
has a critical density of only ~500 cm~^ (Nusbbaumer 



& Storey 1982), two orders of magnitude lower than the 
critical density of [Fe II] and much lower than the electron 
density pertinent to the region of NGC1068 under consid- 
eration here. For these reasons we believe that [Ni II] does 
not affect significantly the measurement of [P II] . 

It is also interesting to note that our data confirm the 



identification of [SIX] by Marconi et al. (1996) and in- 
clude the first detection of [FeXIII] in an extragalactic 
object. The latter identification is also supported by mea- 
surements of the higher ionization green line of [FeXIV] 
(Kraemer & Crenshaw [2000 ). 



3. Discussion 

3.1. [Fell], [PII] and the Fe/P abundance ratio 

The near-IR lines of [P II] and [Fe II] have several interest- 
ing similarities. They lie nearby in wavelength, have simi- 
lar excitation temperatures and critical densities and their 
parent ions have similar ionization potentials and radiative 
recombination coefficients. Using the collision strengths of 
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Table 1. Line fluxes in NGC1068 and template objects 



Line (/xm) 


NGC1068 


Orion 


Relative fluxes'^' 
RCW103 LMC-N63A 


LMC-N49 


[SIII] 0.9529 


1200 


29000 






- 


[CI] 0.985 


72 


20 


— 


— 


— 


[SVIII] 0.9913 


55 


<10 


— 


— 


- 


PaS 1.005 


61: 


1400 






- 


Hell 1.012 


120 


<10 






— 


[SIT] 1.033 


170: 


400 


_ 


_ 


- 


[NI] 1.040 


60: 


14 






- 


[FeXIII] 1.075 


50: 


<10 






— 


Hel 1.083 


1000 


6600 


31 






Pa7 1.094 


110: 


2300 


<15 






[PII] 1.188 


67 


50 


<8 


<6 


3 


Hel 1.197 


<20 


48 


<8 


<6 


<6 


[SIX] 1.252 + Hel 1.253 


72"^ 


66'' 


<8 


<6 


<6 


[Fell] 1.257 


100 


100 


100 


100 


100 


Hel 1.279 


50: 


430 




<6 


<6 


Pa/3 1.282 


230 


4300 


11 


20 


15 


[Fell] 1.321 


33 


34 


30 






[Fell] intensity'^' 


^^20^= 


4.4 


12 


3.0 


4.4 



(1) 



(2) 



Normahzed to I([FeII] 1.257) = 1 00. F luxes for Orion refer to position A of Walmsley et al. (2000), values for the supernova 
remnants are from Oliva et al. (199C) and from the IRSPEC spectra displayed in Fig. |l] The error on line fluxes in 
NGC1068 are tyipically ±10% except for the entries marked with a ":" which are uncertain due to blending. 
Intensity of [Fell] 1.257 in units of 10" 
Contribution from Hel 1.253 should be 
Contribution from [SIX] is negligible. 
The absolute flux calibration is uncertain. 



erg cm s sr 

;10, based on the observed intensities of the other Hel lines. 



Krueger & Czyzak ( |1970| ), Zhang & Pradhan ( |1995D an d 
the transition probabilities of Mendoza & Zeippen ( 1982D , 
Nussbaumer & Storey ( 198^ yields 



n(Fe+) _ 2 -^([FcII] 1.257^m) 
n(P+) ~ ' /([PII] l.mSfim) 



(1) 



which is accurate within a factor of 2 for all the temper- 
atures and densities of interest. The only important dif- 
ference between the two species is that charge exchange 
recombinations between P and neutral hydrogen are 
~2 orders of magnitude less efficient than the Fe++ -I- H° 



reactions (Kingdon & Ferland 1996). This implies 



n(Fe+) n(Fe) 
n(P+) ~MP) 

which combined with Eq. (1) yields 

n(Fe) ^ 2 /([Fell] 1.257Aim) 



/([PII] l.mSfim) 



(2) 



(3) 



i.e. the Fe/P abundance ratio is quite well constrained by 
the [Fell] /[PII] ratio and, if anything, it is overestimated. 
Finally, it is interesting to note that, for a solar Fe/P~100 
abundance ratio, one expects [Fe II]/[P II]~50 i.e. a ratio 
similar to that measured in supernova remnants. 



3.2. Fe/P abundance ratio and dust destruction 

Iron is a well known refractory species whose gas phase 
abundance in the ISM is often found to be down by many 
orders of magnitudes relative its cosmic value. The only 
regions where iron is not found to be significantly depleted 
are those associated with fast (^100 km/s) shocks which 
can effectively destroy the grains by sputtering. In normal 
photoionized regions the depletion of iron ranges between 
the factor of ~0.1 measured in the Orion Bar (Baldwin 
et al. 1996| ) to significantly lower values found in plane- 
tary nebulae (e.g. Oliva et al. 1996, Perinotto et al. 1999| ). 
These differences probably reflect the fact that a variable 
(though small) fraction of iron is locked into relatively soft 
grains which can be easily destroyed without the need of 
fast shocks. 

Phosphorus is a non-refractory species whose mea- 
sured depletion in ionized gas is close to unity. Therefore, 
the Fe/P relative abundance should give a direct estimate 
of the iron depletion or, equivalently, of the presence of 
robust dust in a given region. This is indeed confirmed by 
the data presented here (Fig. |^ and Table |l|) which show 
variations by more than one order of magnitude between 
the large, quasi-solar Fe/P ratio found in SNR's and the 
much smaller values found in the other objects. 
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3.3. [Fell] /[PII] and the origin of [Fell] in galaxies 

As already discussed in the introduction, determining the 
origin of [Fell] Hne emission is of crucial importance for 
any program aiming at using [Fc II] for tracing shock fronts 
and/or constraining the supernovae rate in galaxies. The 
results obtained here indicate that the [Fell] /[PII] ratio 
could provide a clear-cut answer to this problem. The line 
ratio is large (^20) in regions excited by fast shocks while 
low (^2) in normal photoionized region and in NGC1068. 
This indicate that shocks are not the dominant source of 
[Fell] in the central ~200 pc of NGC1068 where, most 
likely, the lines are produced by photoionization from the 
active nucleus, as already indicated by detailed photoion- 
ization modeling (e.g. Kraemer & Crenshaw 2000). The 
copious flux of soft X-rays from the AGN creates a very 
extended partially ionized region which is responsible for 
the strong emission of low ionization species such as [S II] , 
[01], [PII] and [Fell]. However, the latter is relatively 
weak because most of the iron is locked into dust grains. 
The relative intensities of the low and higher ionization 
lines depend on a complex combination of the ionization 
parameter, density and of the spectral shape of the ioniz- 
ing radiation. Nevertheless, the ratio between "close rel- 
atives" such as [Fe 11] / [P II] are little influenced by this 
and, as discussed in Sect. 3.1, almost solely depend on 



the Fe/P relative abundance which, in turn, is a direct 
measurement of iron depletion (Sect. 3.2). 



4. Conclusions 

Given the above results and considerations, we propose 
using the ratio between [Fell] (1.257 /im) and [P II] (1.188 
/im) as a general tool for constraining the origin of [Fe II] 
line emission in galaxies. The diagnostic works as follows 
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- In objects with low [Fell] /[PII] ratios shocks do not 
play an important role in the lines excitation. 

- Large values of [FeII]/[PII] (^20) indicate that the 
emitting gas has recently passed through a fast shock 
which sputtered and destroyed most of the dust grains. 
It is therefore likely that the lines are also produced 
by shock excitation 
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